Realization of a large J2 quasi-2D spin-half Heisenberg system: Li2VOSi04 
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Exchange couplings are calculated for Li2VOSi04 using LDA. While the sum of in-plane couplings 
Ji + J2 ~ 9.5 ± 1.5K and the inter-plane coupling J_l ~ 0.2-0.3K agree with recent experimental 
data, the ratio J2/J1 ~ 12 exceeds the reported value by an order of magnitude. Using geometrical 
considerations, high temperature expansions and perturbative mean field theory, we show that the 
LDA derived exchange constants lead to a remarkably accurate description of the properties of these 
materials including specific heat, susceptibility, Neel temperature and NMR spectra. 
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In many recently discovered magnetic materials the de- 
termination of exchange constants, without input from 
electronic structure calculations, has proven very diffi- 
cult and has often led to wildly incorrect parameter val- 
ues. The interplay of geometry and quantum chemistry 
has yielded many surprises which could not have been 
anticipated without a full calculation. Examples are the 
recently discovered vanadates CaV409 and CaVsOy 
In all these cases the dominant exchange interac- 
tions were resolved and a good understanding of the ma- 
terial properties obtained only after analyses of electronic 
structure calculations were carried out. 

Frustrated square-lattice spin-half Heisenberg antifer- 
romagnets with nearest neighbor exchange Ji and second 
neighbor (diagonal) exchange J2 have received consid- 
erable attention recently. The properties of the model 
with J2=0 (or Ji=0) are well understood at zero and fi- 
nite temperature . The large J2 limit of the model is a 
classic example of quantum order by disorder [^j^] , where 
at the classical level the two sublattices order antiferro- 
magnetically but remain free to rotate with respect to 
each other. This degeneracy is lifted by quantum fluctu- 
ations leading to coUinear magnetic order in a columnar 
pattern. At intermediate J2/J1 there is strong evidence 
for a spin-gap phase, though the nature of this phase is 
not fully resolved yet Q . 

While there has been tremendous theoretical interest 
in these models, there were no known experimental re- 
alizations for intermediate to large J2/J1, until the in- 
vestigation of Li2VOSi04 by Melzi et al. [@,|| Studying 
the splitting patterns of the ^Li NMR spectra, these au- 
thors presented strong evidence for columnar order 0. 
Combining several experiments they derive exchange 
couplings (with J2I Ji ~ 1.1) well into the region where 
model calculations find columnar order. 

However, several puzzling pieces in that excellent and 
detailed study remain: (i) The ratio of exchange con- 
stants was not well determined from the susceptibility 
and specific heat data; we will present electronic struc- 
ture and many-body calculations to show that their esti- 
mate J2/J1 « 1 is off by an order of magnitude, (ii) 



The estimated T — Q moment was anomalously small for 
a system well inside the columnar ordered phase. Taking 
into account the antiferromagnetic inter-plane coupling, 
we propose that the NMR derived moment is small due to 
a cancellation of hyperfine fields from neighboring planes, 
(iii) The order parameter exponent [3 at the transition 
was estimated to be /? « 0.25, which is intermediate be- 
tween 2D Ising and typical 3D exponents. We will show 
that the inter-plane exchange constants differ from the 
largest ones by less than two orders of magnitude. Thus 
a strong crossover between 2D and 3D behavior could 
be expected, (iv) The Neel temperature was nearly field 
independent up to a field of 9T. We will argue that our 
increased estimate of J2 leads to a larger saturation field 
and that combined with non-monotonic dependence of 
Neel temperature on field implies that the experimental 
results are not anomalous. 

Our study of the material Li2 V0Si04 consists of a two 
band tight-binding model fit to the LDA band structure, 
which is then mapped onto a Heisenberg model with in- 
plane (Ji and J2) and inter-plane {Ji_) exchange con- 
stants. Furthermore, we develop high temperature series 
expansions and perturbative mean-field theory for the 
uniform susceptibility and specific heat of the Ji — J2 
model. These allow us to make quantitative comparisons 
with the experiments. 

Li2VOSi04 crystallizes in the tetragonal system, space 
group PA/nmm, containing two formula units per cell 
with a = 6.3682 A and c = 4.449 A. j|] The crystal 
structure of Li2VOSi04 is shown in Fig. ^ The magnet- 
ically active network of spin half ions is built up by 
[V0Si04]^^ layers of VO5 square pyramids sharing cor- 
ners with Si04 tetrahedra, intercalated with Li ions. The 
structure of the V^^ square network suggests, that both 
the nearest neighbor (NN) and the next nearest neighbor 
(NNN) in-plane coupling should be significant, although 
it is at best difficult to decide from general considerations 
which one is dominant. NN coupling is favored by the 
existence of two exchange channels and shorter distance, 
NNN coupling profits from the 'straight' connection be- 
tween pyramids pointing in the same direction. 
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FIG. 1. Perspective view (upper panel) of the crystal 
structure of Li2VOSi04 and projection along [001] (lower 
panel). The VO5 pyramids (large diamonds) share the corners 
of the basal planes with Si04 tetrahedra (small diamonds). 
The Li+ ions are indicated by circles. 

In order to obtain a realistic and reliable hopping part 
of a tight-binding Hamiltonian, band structure calcula- 
tions were performed using the full-potential nonorthog- 
onal local-orbital minimum-basis scheme [ |To| within the 
local density approximation (LDA). In the scalar rela- 
tivistic calculations we used the exchange and correlation 
potential of Perdew and Zunger jl^. V(3s,3p,4s, 4p, 3d), 
0(2s, 2p, 3d), Li(2s, 2p) and Si(3s, 3p, 3d) states, respec- 
tively, were chosen as the basis set. All lower lying states 
were treated as core states. The inclusion of V {3s,3p) 
states in the valence states was necessary to account for 
non-negligible core-core overlaps. The O and Si 3d as well 
as the Li 2p states were taken into account to increase 
the completeness of the basis set. The spatial extension 
of the basis orbitals, controlled by a confining potential 
{r/ro)^, was optimized to minimize the total energy. 

The results of the paramagnetic calculation (see Fig. g) 
show a valence band complex of about 10 eV width with 
two bands crossing the Fermi level. These two bands, due 
to the two V per cell, are well separated by a gap of about 
3 eV from the rest of the valence band complex and show 
mainly V Sd^y and minor 0(2) 2px^y character (oxygens 
of the basal plane of the VO5 pyramid) in the analysis 
of the corresponding orbital-resolved partial densities of 
states (not shown). The valence bands below the gap 
and above the Fermi level have almost pure oxygen and 
vanadium character, respectively. The contribution of Li 



and Si states is negligible in the energy region shown. 




wave vector 



FIG. 2. Band structure and total density of states for 
Li2 VOSi04 (upper panel) and the zoomed bands closest to the 
Fermi level (lower panel). The Fermi level is at zero energy. 
The notation of the symmetry points is as follows: X = (100), 
M = (110), Z = (001) 

The relatively narrow bands at the Fermi level (see 
Fig. ^, lower panel) are half-filled. Therefore, strong cor- 
relation effects can be expected which explain the ex- 
perimentally observed insulating ground state. Because 
the low-lying magnetic excitations involve only those or- 
bitals with unpaired spins corresponding to the half-filled 
bands, we restrict ourselves to a two band tight-binding 
analysis and the discussion of these half-filled bands. 

The dispersion of these bands (see Fig. |^, lower panel) 
has been analyzed in terms of NN transfer ti and NNN 
transfer t2 within the [001] plane (see Fig. |^ lower panel) 
and NN hopping t± between neighboring planes. 

Then, the corresponding dispersion relation of the re- 
lated 2x2 problem takes the form 

E{k) = £0 + 2t2[cos(a;) + cos(?/)] 

±4ti cos(a;/2) cos(y/2) + 2tj_ cos(z) , (1) 

where x = kzO,, y — kyh, z — kzC. 

The assignment of the parameters has been achieved 
by two numerically independent procedures: By straight- 
forward least square fitting of the two bands in all direc- 
tions and by using the energy eigenvalues at different 
selected high symmetry points. The results are shown 
in Table |. The errors can be estimated about 5% for 
the in-plane transfers and 15% for the inter-plane term 
from the differences of both mentioned above fitting pro- 
cedures. These small differences can be ascribed to the 
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influence of higher neighbors. The very good agreement 
of the tight binding fit with the LDA bands justifies a 
posteriori the restriction to NN and NNN couphngs only. 
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0.67 


7.85 
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TABLE I. Transfer integrals of the two-band tight-binding 
model and the corresponding exchange couplings for different 
values of the Hubbard U . 

The resulting transfer integrals enable us to estimate the 
relevant exchange couplings, crucial for the derivation 
and examination of magnetic model Hamiltonians of the 
spin- 1/2 Heiscnberg type: 
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In general, the total exchange J can be divided into an 
antiferromagnetic and a ferromagnetic contribution J — 
jAFM _j_ jFM ^ .j-j-^g strongly correlated limit, valid 
for typical vanadates, the former can be calculated in 
terms of the one-band extended Hubbard model jf^^^ 
= 4:t^/{U — Vi). The index i corresponds to NN and 
NNN, U is the on-site Coulomb repulsion and Vi is the 
inter-site Coulomb interaction. Considering the fact that 
the VO5 pyramids are not directly connected, but via 
Si04 tetrahedra, ferromagnetic contributions are 
expected to be small. For the same reason, the inter- 
site Coulomb interactions Vi should be small compared 
with the on-site repulsion U. From LDA-DMFT(QMC) 
studies and by fitting spectroscopic data to model 
calculations jw), U ~ 4-5 eV is estimated for typical 
vanadates. Therefore, we adopt U =4 eV and U =5 eV as 
representative values to estimate the exchange constants 
and their sensitivity to U. The calculated values for the 
exchange integrals are given in Table |. 

Comparing our calculated exchange couplings with the 
experimental findings we find excellent agreement 
for the smn Ji + J2 = 9.5 ± 1.5 K ||l^ of the in-plane 
couplings, reported from susceptibility data ||] to be 
J1+J2 = 8.2±1 K. In contrast, we find a ratio J2/J1 ^ 12 
which exceeds the experimentally derived ratio in Ref. ||] 
J2/J1 ~ 1.1 ± 0.1 by an order of magnitude. 

In order to understand the experiments better, we turn 
to high temperature expansions for the susceptibility and 
specific heat of the Heisenberg models. Using series ex- 
pansions (r=0) [p^ , non-linear sigma model theory 
(very low-T) , quantum Monte Carlo (QMC) simulations 
(low-T) |17| and high temperature expansions (HTE) 
(high-T) pql, the susceptibility of the nearest-neighbor 
model (Ji = 0) is known accurately for all T. Letting 
J2 = 1 and treating Ji perturbatively, analogous to chain 
mean-field theories, [|l^ leads to the expression 

x(Ji,T)-Xo(T)[l-4JiXo(T)] (3) 



where xo is the susceptibility for the Heisenberg model 
(Ji = 0). As shown in the inset of Fig. |, at T = 
for small J1/J2, this expression compares very well with 
the susceptibility calculated from Ising series expansions. 
|l6| Fig. g also shows that applying Eq. 3 to the finite-T 
QMC data for xo leads to susceptibility values which join 
smoothly with the high-temperature expansion results. 
Thus, we have accurate calculations for the susceptibility 
of the model with small J1/J2 at all T. 
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FIG. 3. Susceptibility (x, with largest x for Ji — 0) for 
J2 = 9K, g = 2 and J1/J2 = 0, 0.1, 0.2. The low temperature 
data is obtained from QMC combined with Eq. 3, while the 
high temperature data comes from HTE. The inset shows 
Ising series expansion calculations and Eq. 3 for T = 0, J2 = 1. 

Rather than find a fit for the exchange constants, in 
Fig. H we show the susceptibility with g = 2; J2 = 9K; 
and J1/J2 = 0,0.1 and 0.2. The results are close to 
experimental values Q . We note that the agreement will 
be improved by going to the lower limit of the calculated 
exchange constants and slightly larger g- values. 

The specific heat data was the primary source for the 
J1/J2 ~ 1 conclusion by Melzi et al. g. They found 
that the peak value of the specific heat in Li2VOSi04 
was 0.436(4)i? at T™ = 3.5(1). We find that for the 
pure Heisenberg model the specific heat peaks at T™ — 
0.60(4) J with a peak value of 0.455(10)i?, in agreement 
with Ref. [^. With small J1/J2 the peak shifts to lower 
temperature and the specific heat becomes flatter. The 
fact that the values for the pure Heisenberg model are 
close to the experiments strongly favors a small Ji . 

One of the most puzzling aspects of the experimental 
results is the small moment of 0.24 /ib at T = 0, 
obtained from the NMR split patterns. In contrast, the 
moment of the square-lattice Heisenberg model is well 
known to be ~ 0.6fiB- Taking into account the con- 
siderable antiferromagnetic inter-plane coupling Jj^ re- 
sulting from our calculation, a part of the discrepancy 
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could be understood: The Li nuclei sit off-center be- 
tween two planes, which results in a partial cancellation 
of the hyperfine fields fi'om antiferromagnetically ordered 
NN and NNN V sites in neighboring planes (see Fig. |^). 
This partial cancellation does not change the arguments 
of Melzi et al. for the pattern of line-splitting (including 
intensities) and its relation to columnar order because 
the ordering pattern inside the planes remains the same. 
However, it leads to a reduction in the effective hyper- 
fine coupling and hence to an enhancement of magnetic 
moment derived from the line shift. Taking into account 
the calculated two center overlap integrals for Li and NN 
and NNN V 3d orbitals, respectively, (see Fig. ^ a crude 
estimate from Slater-Koster integrals suggests that the 
NMR split would be reduced by an additional factor of 
about 2. This results in a moment of about 0.5 /x^ much 
closer to the value expected for the 2D Heisenberg model. 



change constants for the material Li2VOSi04 and devel- 
oped numerical studies for the Heisenberg model to show 
remarkable consistency with many experimental proper- 
ties. Electronic structure calculations on the closely re- 
lated material Li2VOGe04 will be presented in a forth- 
coming publication. The key differences are the consid- 
erably smaller J2/J1 ratio and coupling to higher neigh- 
bors in Li2VOGe04. Finally, we note that both these 
materials have a substantial 3D coupling, which leads to 
long-range order at finite T. It would be interesting to 
find a material with large- J2 that was nearly 2D, thus 
closer to exhibiting purely quantum order by disorder. 
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Deutsche Forschungsgemeinschaft, and by the Australian 
Research Council. 
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FIG. 4. Sketch of the different magnetic environments for 
the ^Li NMR. The Li and V sites are represented by black 
and gray circles, respectively. The arrows indicate the direc- 
tion of the V spin. Full lines symbolize the stronger inter- 
action with the upper plane (NN), dashed lines the weaker 
interaction with the lower plane (NNN). The different envi- 
ronments cause (a) no NMR shift due to complete in-plane 
moment-cancellation; (b) up or down shift with partial in- 
ter-plane moment-compensation. 

We now turn to the inter-plane couplings and the mea- 
surements of the Neel temperature, Tn- Applying the 
expression Tn « 0.36 Jj_^^(T/v) |^ (^ is the in-plane cor- 
relation length), to our calculated exchange constants, 
leads to the estimate T/v « 3.6 ± 0.4 K, which is remark- 
ably close to the experimental value of 2.8 K. Further- 
more, the saturation field for our calculated exchange 
constants is about 30 T, which is much bigger than the 
9 T field applied by Melzi et al. The Neel temperature 
should go to zero at the saturation field. However, we 
note that due to suppression of spin fiuctuation the Neel 
temperature can increase slightly with field, as happens 
in the purely 2D model. Thus, the experimental result of 
very weak field dependence of the Neel temperature up 
to 9 Tesla is consistent with our expectations. The ap- 
preciable but still small 3D couplings should also give rise 
to 3D critical behavior at the finite temperature transi- 
tion with strong crossover effects. These results on the 
field dependence of the Neel temperature and the critical 
behavior at the transition in weakly coupled Heisenberg 
systems deserve further theoretical attention. 

To summarize, we have used LDA to calculate ex- 
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